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Abstract: In this paper we study the reductions of evolutionary PDEs on the manifold of the stationary 
points of time-dependent symmetries. In particular we describe how that the finite dimensional Hamiltonian 
structure of the reduced system is obtained from the Hamiltonian structure of the initial PDE and we 
construct the time-dependent Hamiltonian function. We also present a very general Lagrangian formulation 

^ ' of the procedure of reduction. As an application we consider the case of the Painleve equations PI, PII, PHI, 

PVI and also certain higher order systems appeared in the theory of Frobcnius manifolds. 
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Introduction 

Consider an evolutionary PDE with one spatial variable 

ut^F{u,u^,...,u^"'^) (0.1) 

and a symmetry of this equation, i.e. another system 

u, = GKm„...,uW) (0.2) 

commuting with the first one, 

The set of the stationary points u^ = of the symmetry is a finite-dimensional invariant manifold for 
the system (0.1). Particularly, in important examples, the invariant manifold can be described as a set of 
stationary points of a first integral of the system (0.1): 

" 0, 



Su{x) 

I^ f L{u,u,,...,u("^)dx, ^=0. 

In this case it is known ([BN], [Mo]) that the restriction of the initial PDE to the invariant submanifold is 
a Hamiltonian system of ODEs. In particular Bogoyavlenskii and Novikov [BN] found a universal scheme to 
construct the Hamiltonian function of the reduced system in terms of the Hamiltonian of the original PDE. 
In this paper we extend this scheme to more general finite dimensional invariant submanifolds specified 
by local x- and time-dependent symmetries and conservative quantities of the evolutionary equation. To 
distinguish this class of symmetries from the previous one we will call them scaling symmetries. We show 
that the restriction of the starting equation on the finite dimensional manifold admits a natural description 
as a Hamiltonian system with time-dependent Hamiltonian. 

The best known class of examples of evolutionary PDEs admitting nontrivial symmetries and 
conservation laws are integrable systems of soliton theory (sec [SM] and references therein). The finite 
dimensional manifolds of the stationary points of integrable systems are typically described by ODEs of 
Painleve type [AS] , [CD] . For the simplest examples of these restrictions the Hamiltonian structure is already 
known. For example for the classical six Painleve equations the Hamiltonian description was found by 
Okamoto, [O]. Although the relationship between the starting PDE and the reduced ODE is clear and 
has been investigated quite a lot (see, e.g. [AC], [AS]), the relationship between the starting Hamiltonian 
structure and the reduced one has not been elucidated. This work will give a contribution in understanding 
of this relationship. 

As a first result (see section 2) we prove that the finite dimensional Hamiltonian structure of the 
ODEs is obtained from the Hamiltonian structure of the starting PDE, via scaling reduction. Particularly, 
we construct the time-dependent Hamiltonian function of the reduced system. In the time-independent 
case this procedure coincides with the well known stationary-flow reduction discovered by Bogoyavlenskii 
and Novikov [BN]. As an application we present the case of PI, PII, PHI, PVI and also certain higher order 
systems appeared recently in the theory of Frobenius manifolds [Dl]. 

As a second result (see section 3) we present a very general Lagrangian formulation of the procedure 
of reduction of an evolutionary system (0.1). Namely, we prove that this restriction is again a Lagrangian 
system with the Lagrangian function A, such that 

dA _ dL 
dx dt 

The work is structured as follows: after recalling, in Section 1, some basic facts about the Hamiltonian 
structure of the evolutionary PDEs, and briefly summarizing the method of reduction of evolutionary flows 



on the manifold of stationary points of their integral, introduced by Bogoyavlenskii and Novikov [BN], in 
Section 2 we consider the generalization of this procedure to scaling symmetries. The reduced flow is a time- 
dependent Hamiltonian system, and in Theorem 2.1 we give the relationship between the infinite-dimensional 
Hamiltonian structure and the reduced one. 

Section 3 is devoted to a Lagrangian approach to the problem: after describing the general framework, 
in Theorem 3.1 we give the procedure of reduction and we construct the reduced Lagrangian function. In 
Section 3.2 we establish the relationship with the Hamiltonian approach. As an application we study the 
Lagrangian reduction of KdV on the manifold of the fixed points of the 7-th flow. 

Section 4 contains the application of the theory to the scaling reductions from KdV, niKdV and Sine- 
Gordon equations respectively to Painleve I, Painleve II and III. These examples are studied both from the 
Hamiltonian and the Lagrangian point of view. 

In Section 5 we study the n-waves equation and his scaling reduction to a system of commuting 
Hamiltonian flows on the Lie algebra so(n). The reduced system is a no n- autonomous Hamiltonian system 
w.r.t. the Poisson structure of so(n). In particular, for n = 3, adding an additional symmetry condition, 
one arrives at Painleve VI equation. 



1. Infinite dimensional Hamiltonian structures and stationary flows reduction 

Let us consider the phase space dJl of smooth maps of the circle into some smooth n-dimensional manifold. 
Actually we can forget about the boundary conditions when dealing with local functionals only. We denote 
by ^ the space of smooth functionals on 93t of the form 

(m)( 



F{u) ^ / f{x,u{x),u^{x), . . . ,u'^"^>{x))dx, 

where the density / depends only on a finite number of derivatives of u. On the space ^ the variational 
derivative g^i(^\ is defined by 

f 5F 

6F^ / ——-du'^(x)dx. 

Explicitly, 

^F df V-, -,,'» d' df 



Y.(-^) 



Su^{x) du'^{x) ^-^ dx'° du'^^\x) 

One can define on Wl the (formal) Poisson brackets 



N 



{u\x),u^{y)} ^w\x,y) ^^X^S {x - y) 



fe=o 



where A^. depends on a finite number of derivatives of u. This induces on ^ the Poisson bracket 



where 






dx 

k=0 



A Hamiltonian system on $H has then the form 

s:tt 

ulix)^iu^ix),H}^P^^-^^. 



In particular, we consider so called Gardncr-Zakharov-Faddev bracket P^ ~ '^*"'^- ^^ ^^^^ '^^^^ ^ 
Hamiltonian system has the form 

..(.) = M.),i/} = A_^ (1.1) 



with Poisson bracket 



Let us consider a first integral 



\FG}= f ^^ — -^dx 
J Su(x) dx du(x) 



I = L{x, u{x), Ux{x), . . . , u^"^)dx, 
where L does not depend on t. The generalized Euler-Lagrange equation 

SI 



6u(x) 



(1.2) 



generically is a ODE of order 2n fixing the 2n-dimensional manifold S of the stationary points of the first 
integral /. Because of the Lax lemma (see [Mo]) this submanifold is invariant w.r.t the evolutionary equation 
(1.1). The functional L is the Lagrangian of the x-flow defined by (1.2). If L is nondegenerate, then it defines 
also on S the natural system of canonical coordinates 

q^^u'^'-'^\ i^l,2,....,n 
SI 

P^ 



and equation (1.2) can be put in the Hamiltonian form 

where H is the generalized Legendre transform of L: 

SI 



h=-l+j: 



Sui"- 
1 



u« 



which, in terms of the canonical coordinates takes the form: 

" — P'ii- 

It is well known that the starting PDE can be restricted on 6 and the restriction is a Hamiltonian system 
of ODEs. In particular Bogoyavlenskii and Novikov discovered the algorithm to construct the Hamiltonian 
functions of the reductions in terms of the Hamiltonian of the original evolutionary equation. They considered 
the case of a hierarchy of evolutionary equations 

du d SI, 



dt,. dx Su{x) 

with If, = J Li_ {u{x), Ux{x), . . . , u^"'''{x))dx, and they described the reduction procedure of the fc-th flow on 
the finite dimensional manifold of the stationary points of the j-th fiow. They proved that all the flows of 
the hierarchy reduce to finite dimensional Hamiltonian system. The Hamiltonian function for the reduced 
k-th fiow, {—Qk.j), is determined by : 



5L d SL d 



Su{x) dx Su{x) dx 
Mokhov [Mo] generalized this result to not necessarily Hamiltonian evolutionary PDEs. 



2. Scaling reductions of evolutionary systems: 
Hamiltonian formulation 



In this Section we extend the Bogojavlenskii-Novikov scheme to finite dimensional invariant submanifolds 
specified by time-dependent local symmetries. 

We start from a partial differential equation of order m on the functional space 971, describing the 
evolution of the function u{x) in the time t and a scaling symmetry 

Us = G{x,t,u,Ux, ■ ■ ■ ,u''''^)- 

Our main assumption is that the set of stationary points of the symmetry can be formally represented in 
the Euler-Lagrange form 

1= L{x,t,u,Ux,...,u^"-'>)dx, ~u—^- 

It is an ordinary differential equation of order 2n depending explicitly on the parameter t. If L is 
nondegenerate, the space of the solutions is a 2n dimensional manifold &, which naturally carries a system of 
canonical coordinates. As in Section 1 we will show that, in these coordinates, the Euler-Lagrange equation 
is Hamiltonian, with Hamiltonian function H , obtained from L via Legendre transform: 

Following the scheme of [BN], we prove that one can reduce on 6 also the equation of the evolution in i, 
which results to be a Hamiltonian system. We also give a universal scheme to produce the time-dependent 
Hamiltonian function of this reduced system. Indeed the following theorem holds: 

Theorem 2.1: Ij the evolutionary PDE: 

Ut^F{u,u^,...,u^"''>), 

admits a nondegenerate scaling symmetry, then, on the manifold & of the stationary points of the symmetry: 

SI 

0, 



Su{x) 

/T T 
L(x,t,u,Ux, ■ ■ ■ ,u )dx, — = 0, 

dt 

it reduces to a Hamiltonian motion in t, for the time dependent Hamiltonian function {—Q), that is the 
reduction on 6 of 

Q-^-J2p^^^ (2-1) 

1=1 



(2n~ 1) 

where pi, qi are the canonical coordinates on &, expressed in terms of u, u^, . . . ,u , and the function A 

is determined by 

^ = ^. (2.2) 

dt dx ^ ' 

Proof: We prove the theorem in three steps: first we describe the submanifold S of stationary points of 
the symmetry /, where we introduce a system of canonical coordinates; then we deduce, on 6, a zero- 
curvature equation for (— Q) and the Hamihonian function H of the reduced a:;-flow. FinaUy we prove that 
the restricted i-flow is Hamiltonian on S,with Hamihonian function (— Q). 

1) The manifold 6 is the 2n~dimensional manifold of the solutions of the Euler-Lagrange equation 

51 



= 0. (2.3) 



5u{x) 
It is invariant under the t-flow and it naturally carries a system ofcanonical coordinates: 

ci^=u'''-^\ i^l,2,...,n (2.4a) 

obtained via generalized Lagrange transform (here we suppose that the generalized Lagrangian L is 
nondegenerate). Observe that now the pi depend on x and on t. 

Reversing relations (2.4), one can express the derivatives u,Ux, ■ ■ ■ ,u^'^"^^' in terms of the canonical 
coordinates pi and qi, x and t; explicitely: 

w'"-' == (qJx = gi{x,t,q^,...,qn,Pn) 

y_{n+l) = g^(^x,t,q^,. .. ,qn,Pn,Pn-l) 



y(2n-l) ^ g^(^^^ t,q^,..., qn,Pn, • • • ,Pi)- 

Observe that (2.4) gives the identities: 

yPiJx + gq^ — Su 
dH 



(pX + ^^O, i>l (2.5) 



where H is the generalized Legendre transform of L: 



^I (^) 



2 — 1 



which, in terms of the canonical coordinates takes the form 

n 



dq- 
dx 



The first of identities (2.5) allows us to express the higher derivatives u^™^ for to > 2n in terms of x,t, 
Pi, qi and p['^ with / = 1,...,to — 2n+l, explicitly: 



u 



^^"^ = 9n+l{x,t, q,,..., qn,Pn, ■ ■ ■ ,P^,{Pi)x) 



U ^ gm-n+l[X,t, q^,...,qn,Pn, ■■■,?!, ■■■,[Pl) )■ 

6 



On & it reduces to {Pi)x + § — = 0, and the system (2.5) is a canonical Hamiltonian system, with Haniiltonian 
function H, giving the reduced x-flow. 

Now we will show that also the i-flow reduces on 6 with Hamiltonian function (— Q). 

Firstly we observe that Q is a function of x,t, u and its x-derivatives up to the order {m + n), then 
it can be rewritten in terms of x, t, {pi, qi) and pi'' up to the order I = m — n + 1. 

We denote with / a function f{x, t, u{x), . . . , u{x)'-^^) reduced on 6; notice that, if j > 2n, then the 
reduction can be done using the relation 

Then / does depend explicitly only on the pi and qi, for i ~ 1, . . . , n and on the time t. In fact differentiating 
(t.6) one obtains the derivatives pi'' in terms of the canonical coordinates {pi,qi). 

2) We consider the derivative 

dL dL s—^ dL dqi ■^-^ dL dpi 

'dt ^~dt ^ ^^dq^'dt ^ ^^djH~dF ^ 

_ _ dH_ ^^ dH_d(ii_ .^ d^q, 

" dt ^ dq, dt^^^'dxdf ^ ' 

i=l ^* i=l 

From the fact that / is a first integral, one deduces that ^ must be the total derivative in x of a functional 
A that does depend on x,t,(pi, qi) and pi'' up to the order ^ == tti — n + 1; we have: 

dA _dQ -^ dQ dqi J^ dQ dpi "^A dQ d (^j 
dx~ dx^^ % dx^^ dp, dx ^^ d(p,Y') dx^^'' 

n , , n ,2 

>-^ dp, dqi -s-^ d qi 



dx dt ^-^ dxdt 

i—l i—1 

_dQ ^ dQdH_ _^ dQdH_ ""^^ dQ d (,) 
" 9x ^ -^ % dp, y^ dp, dq, ^ -^ 9(pJ(*) dx^^'' ^ 

-A dH dqi A^ d2 dqi dQ 

-Ea^^ + E^^^ + ^(^.)^ + a^(p.)^- (2-8) 

Then equation (2.2) gives: 

dH dQ ^dQdH_ ^ dQdH_ dQ 
dt dx ^ dqi dpi -^ dpi dqi dpi ' 

h d{p,Y^) dx^^^' + dt^^''^^ dqi dt ~^' 

which can be rewritten as 

dH dqi ( , , dH\ d ^ ,„ „, 

At this point we need the 

Lemma 2.1: On the submanifold & the following relation holds: 

7 



Proof: See Appendix 2. A 

Hence, on the submanifold & cq. (2.9) reduces to: 

dH dQ ^ dQ dH ^ dQ dH 
dt dx ^ dqi dpi -^ dpi dqi 

This is a zero-curvature equation: 

This completes the second step in the proof of the theorem. 

3. Now we win construct the Hamihonian system inductively; to this end we need a further lemma: 
Lemma 2.2: The fundamental relation 

dt ~ dp, ■ ^ > 

holds. 

Proof: See Appendix 2. A 

For simplicity, here and in the following we omit the "tilde" sign: Q will indicate the reduced function 
on &. Now, we assume that -^ — —g^ — ^{qi,Q} and we prove inductively that the same relation holds 

for q-_^-,. The scheme of the procedure is the same as in [BN], the only differences are the contributions of 
the partial derivatives in t and x. Indeed, 

Using the Jacobi identity and the zero-curvature relation we get 

—T^ = ^{"^'5'.} - {<li+i^Q}^ i = 1,2, ...,n- 1. 

Here the term {^,(?j} is zero for every i ^ n since {%;) — ^^- Indeed L depends on u and on the 
derivatives of u up to the order n. This means that, restricted on 6, it depends on g^ , g^ , . . . , qn+i- Then, 
there is no dependence on the p. for i ^ n. Finally we get 

dH ,/dL\ 

|-^,ft}=0 i<n. 



Hence we have proved that 



-rf = -{9,,^}, 2 = l,2,....,n. 

at 



Now we prove that -J^ = ~{Pn Q} by induction, starting from p„. 

8 



This comes from the eomm.utativity of the flows 



d d d d 

dt dx dx dt 



exphcitly: 



d_f d_ \ _ d rOH 
dt \ dx "/ dt \ dpn 

_d dH ^ dH dQ d^H d 
dt dpn ^ dpndq. dp^ dpi dt 

dH ^ dH dQ d^H d 
dt ^ dpndq. dp. dpi dt 

On the other hand, using the Jacobi identity and the zero-curvature equation, one can write 

d fdqn\ d dQ 

do 
= - {qn, {Q, H}} - {{qn, H}Q} - {qn, -^} - 

dH dH 

dH y. dH dQ d^H dQ 
dt ^^dpndq^dp. dp^ dq„ 

Comparing the two expressions and noticing that S-^ ^ because of the nondegeneracy, we get 

dPn r ^. dQ 

- Now we suppose that —rf- — {p^, Q} and we deduce the same for p._i. Indeed, 

±(dpA rdQ ^.^d_dQ^ 
dx\dt) dq. ' dx dq. 

= -{fe,0},H}-fe,g} = 

fifj 

= -{fe,i/},Q} + fe,— } = 

"-1 02 TT an pa 



and 





, on , v^ 


d (dp\ 
dt\dx ) 


d (dH\ 
dt \ dq. J 

d dH ^ 

f)f dn ^ 



d'H dQ d^H dQ d'H dQ 



dp^ dq.dpndqn dq^dp^_^ dq^_^ 



d^H dq^ d^H dpn d^H dp^_^ 



9tk '^1^ ~{ dq.dq. dt^ dq.dpn dt^ dq^dp^_^ 

"^1 a2TT ari Pfiu ^/O Fii} 



dt.^ 



dH ;-A d^H dQ d'H dQ d'H dp^_, 

"' dt^ ^^dq^dq^ dp^ dq.dpndqn dq,dp^_-^ dt^ 

9 



where g-^ = 1 • Comparing the two expressions we get — ^^ — g ^ ; hence it follows that 

^-^-fpQ} z-1 2 n 



Q.E.D. 



Remark: The definition of Q: 

2—1 

looks very similar to the definition of the Hamiltonian function H of the x-flow: 



i—l 

Here a symmetry between x and t seems to appear: one could be tempted to read the definition of Q as 
a Legendre transform and hence to read A as the Lagrangian of the t-flow. But it is not completely true: 
indeed the coordinates qi and pi are obtained from the Lagrangian L, they are not, a priori, good coordinates 
for A. In the next chapter we will perform a change of coordinates on 6, in order to read A as Lagrangian 
function. 

2. A Appendix 

Proof of Lemma 2.1: We observe that the recursive relation 

dQ \ d f±r)\_±( dQ 



,9(pJ(j-i)y d{p^Yj)\dx^J dx\d{p,)(i') 
holds for j > 1. Indeed 



d{p,)U)\dx^j ^v%5(pjw;^^*^" ^^\dp^d{p,Y^) 



When we reduce on &: 

dQ \ d f dQ 



.d{p,Y^-^'>J dx \d{p,)(-i') 

But Q depends on (Pip-'' up to a finite order, then 

^ ' V,>1. 



a(pjo) 



Q.E.D. 



Proof of Lemma 2.2: The expansion of -^Q in powers of (Pi)x near the point ^^ reads 



d 
dx 



d{pjx \dx 



, , dH\ f ^ ^ dH 



2 



10 



The zero order term is 



d 
dx 



dH 



by virtue of the zero-curvature equation. The first order coefficient is 



dipjx \dx 



= 


\d( \1 


\dpj 

dp, U UbJf') dpi 


dx\d{p^)x'^ ) 




dx\d{p^)x'^ ) 



where the only non zero term is g^, by virtue of Lemma 2.1. Hence we obtain the power series 
expansion oi j-Q up to the first order: 

dH dQ f , , dH 



Q.E.D. 



this, compared with the left-hand side of eq. (2.9), gives the relation (2.f f ). 



3. Scaling reductions of evolutionary systems: 
Lagrangian formulation 



3.1 General framework 

The basic idea is to develop a reduction method dealing on the same footing with x and t. The starting 
point is always the evolutionary PDE 



ut^F{u,Ux, ,«(")) 



(3.1) 



in the space 9Jl described in Section f .1. The first step of our construction is to read u as a function of x 

and t and to consider equation (3.1) as a definition of ^('"^(a;, t) in terms oi u{x,t),Ux{x,t), ,u^"^^^\x,t) 

and ut(x,t). 

This corresponds to consider as "coordinates" in 971, instead oi%i{x,t) and its derivatives in x: 



lA, llx 1 Uxx^ • • • 



(here and in the following u indicate the function u{x,t)), the system 



%Ux, 



,U^ ',Ut,Uxt,-- ■,ul ',Uu,- 



By virtue of the reversibility of (3.1) in u''"^'>{x,t) it is possible to perform this "change of variables" . If one 
introduce the vector 

the new system of "coordinates" in 9Jl is given by u{x, t) and its derivatives in t: 

U,Ut,Utt, ■ ■■ 

11 



At this point one takes a first integral of eq (3.1), i.e. a functional 

/= /i(a;,i,u(a;),u^(x),...,u(")(x) jdx (3.2) 

in the space $H, such that 

-^ = 0. (3.3) 

du[x) 

This Euler-Lagrange equation defines a finite dimensional manifold S, i.e. the set of the fixed points of 
/. Indeed the Euler-Lagrange equation (3.3) is an ODE of order 2n, so that the space of the solutions is 
a 2ri-dimensional manifold; 6 is modeled on this space, having as coordinates certain combinations of the 
initial values, i.e. of the first (2n — 1) x-derivatives of u{x) evaluated at xq- 

In Lemma 3 . 2 below, we rewrite the definition of the manifold 6 in terms oiu{t),ut(t), . . . and of a functional 

J= / A(x,t,u{t),ut{t),...,u'''^\t)]dt, 

where A can be calculated from L (see eq. (3.4)), and the order (3 of derivation in t depends on the ratio 
between m and n, as we will show in detail in Section 3.3. 

In Theorem 3.1 we will prove that A I x, t,u{t), ut{t), . . . ,u^^^\t) j is the generalized Lagrangian for the i-flow 

reduced on 6. Indeed equation (3.1) can be rewritten in form of a Euler-Lagrange equation: 

5J 



Suit) 
for the vector 



= 0, 



6J ( 5J 6J SJ 



5u{t) \Su{t) ' 5u^{ty"' ,5m(™-i) {t) 

where 

5J _ dA ^ c d" dA 

i 
Explicitly, equation (3.1) reads 



In the multiindex (i, a) the Latin character indicates the order in the x-derivative, the Greek indicate the 
order in the i-derivative. 



dA ^ y. d" dA 
dA ^ . d" dA 

du,+2^^ > dr d^i-)-^ 



d A ^ ^ rf° dA 



au(m-l) A^"- > df Q^(m-l,a) 



We will formalize these facts in the following Theorem 3.1; here we give an idea of the proof, ignoring all the 
calculations, that we will concentrate in Lemma 3.1 and 3.2, the proof of which is postponed in Appendix 
3.A. 

The proof is performed in four steps: firstly we define the new system of "coordinates" in 9Jl, and 
we give some useful relations between the new and the old "coordinates" . As a second step we rewrite the 

12 



Lagrangian density L[u(x)j in terms of the new "coordinates" and we construct, starting from L[u{t)j, the 
Lagrangian density A(u(i)). 

The third step consists in recovering the relation between performing the variation of L (in x) and of 
A (in t). The most relevant relation is that given in Lemma 3.1. This relation is necessary to rewrite the 
Euler-Lagrange equation (3.2), defining 6, as a condition on A. The explicit form of this condition is given 
in Lemma 3.2. 

Finally we prove that, under this condition, i.e. after performing the reduction on 6, the starting 
evolution equation (3.1), reads as an Euler-Lagrange equation for A. 

The method of Hamiltonian reduction described in Chapter 2 allows us to put a canonical system of 
coordinates {pi, qi} on 6 (see formula (2.4)). These coordinates are obtained from L via generalized Lagrange 
transform, so that they are, in a certain sense, adapted to the a;-flow. This means that in these coordinates 
the reduced x-flow is a Hamiltonian system. Theorem 2.1 also gives the explicit form of the Hamiltonian 
function 

H = -L + y^^pi{qi)x. 

i 

The method of Lagrangian reduction which we describe in this Chapter, still allows us to define a system of 
canonical coordinates: we will call it {pi,qi}. These coordinates are obtained from A, i.e. they are adapted 
to the i-flow; in fact we will prove (see Section 3.3) that, in these coordinates, the reduced i-flow is a 
Hamiltonian system, with Hamiltonian function 

-Q^-k + ^p,{q,)t. 

i 

When rewritten in terms of {pi,qi\, the Hamiltonian Q coincides with the Hamiltonian function Q 
constructed by Bogoyavlenskii and Novikov. 

In this sense the alternative definition of Q given by us in Theorem 2.1: 

-Q = -K + ^p^{qi)t, 

i 

is a Legendre transformation, if one uses the right system of canonical coordinates (see below) . 

3.2 Lagrangian reduction 
Theorem 3.1: If the evolutionary PDE: 

ut = F{u,u,, ,u(™)), (3.3) 

admits a nondegenerate scaling symmetry, then, on the manifold & of the stationary points of the symmetry: 

61 



Su{x) 



0, 



/T T 
L{x,t,u,Ux,. ..,u )dx, 'Jf = ^' 

it reduces to a Lagrangian motion in t, for the time dependent Lagrangian function A, determined by: 

dt dx 
13 



Proof: (in the following we consider the case m < n < 2m. The same holds in the case (a ~ l)m < n < am, 
as we will show in Section 3.3). We prove the theorem in four steps: 

1. Change of "coordinates": Let us assume that the evolutionary equation (3.3) depends on u{x) and on 
its x-derivatives up to finite order m,, and that this equation is invertible in w'™'. In this case we can read 
(3.3) as a definition of u*^™-* in terms of u, u^, . . . , w'^™^^^ and Ut- 



i^"'^^.fo(u,u,,...,u^"'-^\ 



Ut 



(3.5) 



Differentiating cq. (3.5) in x one obtains all the a;- derivatives of u of order greater then m in terms 
oi u,Ux, ■ ■ ■ , u^"^^^' and their t-derivatives: 



,Ut,Uxt ••■,"( 
(m-1) 



(m-l)N 



, m(™+") = fn{u,Ux,. . . ,U^'"'-~^\ut,Uxt ■ ■ ■ ,u]^'^\uu,UxU, ■ ■ ■ ,Ult 

Explicitly, the first relation has the form 



(m+l) _ "" I '^" 






9a; 9u Su^^^^) ' i9ut 

and in general, using the multiindex notation introduced in Section 3.1, 



-Wxt, 



9a; 



au(fe,/3) 



fc=0 /3=0 

where the higher order a in the i-derivative is fixed by [a — l)m < j < am. 
This completes the construction of the map from 



,1^ , , 



to the new system of "coordinates" 



U,U^,...,U^ >,Ut,...,ul ',Utt,Uoott,-- ■,ult 



Here below we list some noteworthy relationships between the two system (they will be useful in the 
following): 



dt 



du 



-Ut 



Ouim-l) "* 



dut 



-Utt, 



(3.6a) 



dx 
d 



dut 



d /du(''> 






,(fe) 



9^.^'=-^) 



9uJ >^"t 









au 



(3.66) 
(3.6c) 

(3.6d) 
(3.6e) 
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2. Lagrangian densities: The Lagrangian L defining the symmetry, depends on u,Ux,. ■ ■ jw'"-*, so that its 
derivative ^ depends on u,Ux, ■■■ , u(™+") . In terms of the new "coordinates" one may rewrite L as 

L(x,i,«,u,,..., «("-!), ut,..., 7/1"-"') (3.7) 

and 

+ ( -Q-(^ j Wi"^ (u, . . . , Mtt, . . . , M^r^"') (3.8a). 

Of course, (3.8a) coincides with 

dL dL dL dL fm-i'i dL dL in-m\ ,^^,^ 

-17 = ^ + -Tr^t + ■■■ + -m — TT^* + -^^^tt + ■■■ + —, tm ', 3.86 

dt dt du 9u(™-i) * dut du^"''"'^ tt ^ 



where 



dL fdL\ ^ f dL \du^''^ 



dt \dt I ^^ \ dui^) ) dt 

k—7n 



dl f dL\ ^ ( dL\ du^^) • n 1 ex oh^ 

Z. brm h;7iT ^-0,...,m-l (3.96) 



k—m 



dL ^ f dL \du^''^ 

E hro h-TI) ^ = 0,...,m-l. (3.9c) 



From the fact that / is a first integral, it follows that there exists a functional 

A(x,t,u,Ux,...,u'- ',ut,...,ul ',uu,---,ult '), 



such that 



where 



dL _ dk 
dt dx 



dA_dK dk dk („_i) 

dx dx du 9u(™^^) 

dk f™-,, crn-n x dk dk (m-i) 



|9A („) \ , dk dk („_„) 



H ^^tt 



3. Variations: Our aim is to reduce equation (3.3) on the space 6 of the stationary points of / = J Ldx. 
This finite-dimensional manifold is defined by the Euler-Lagrange equation 

'' 0. 



5u{x) 
15 



This is a variational equation in the old "coordinates" u, w^;, • • •; how can we define the same manifold 6 in 

SI 

7{xj 



terms of the new "coordinates"? We must express gMpr in terms of A and its variations. To this end we 



recall that 

SI _ dL ^ ■ dJ dL 

and wc first express the terms 

5L 



for J > in terms of A and the the new "coordinates" , namely: 
Lemma 3.1: The following recurrence relation holds: 

dA d ( dk\ _( TT \ ^ ( IT \ du(^-'^ . 

Proof: see Appendix 3. A 

The proof of Lemma 3.1 is based on the comparison of (3.8) and (3.10) and their partial derivatives w.r.t. 
u/ and ull ■ With a similar technique, and using equation (3.11), one can proves the fundamental 

Lemma 3.2: The (generalized) Euler-Lagrange equation 

SI 



Su{x) 
is equivalent to the condition: 







dA d f dA , , , 

= 0. (3.12) 



Qu{m-1) diV5u(™-l) 



Proof: see Appendix 3. A 
Introducing the functional 

J = J A(^x,t,u{t),u4t), . . .,u^"^-'\t),ut{t), . . . ,uir'""'^(t) ]dt, 

equation (3.12) reads 

'' 0. 



(5M("-i)(f) 
Notice that the object in the left hand side is the last component of the vector 

SJ ( SJ SJ SJ 



5u{t) \Su{t) ' 5u^{ty'' ,5m(™-i) {t) 



SJ 
ncic wc prove luai an iiic couipuncnis ui uiie vecuui ■ 

can be done recursively, by mean of 



4. Reduced evolutionary equation: Here we prove that all the components of the vector ^-i ■. are zero. This 



16 



Lemma 3.3: The following recurrence relation holds: 
SJ 6J 9m("') < 



(5u(*-i)(t) (5M(™-i)(i) 5uW dxySui^^it) 



SJ 



i = l, 



(3.13a) 



Proof: see Appendix 3. A 

Indeed, Lemma 3.2 states that the {m — l)-th component ^ (m-i)(f\ is zero when reduced on &, hence, by 
virtue of (3.13a), all the components of ^i^Ls vanish on 6. 
This is the Euler-Lagrange equation for the Lagrangian 

A(x,t,u,Ux,...,u ,ut,...,ul ',utt,---,ul^ '). 

Q.E.D. 



Remark: Equation (3.13a) can be rewritten as 



_di 


dA d r dk \ d^ 

.(-1) dt[d^f-^) J ^ df- 

' dk d f dk \ 


( ^^ V 

'\duf-'^)\ 
du^"''> d 
du^^^ dx 


■ dk 

dui'i 


d 
" It 


^ dk ^ 

Kdu^^^j 


d^ 
+ dt^ 


( dk 




[a«(™-i) dt[ou("^-^))\ 


yd^ 



for i=l,...,?7i — 1. And The Eulcr Lagrange equation reads 



dA 



d dA 



(9u(^) dt Q^( 

dA _ d I dA \ , d^ dA _ n. 



= i = n — m, ...,m— 1 

i = 0, . . . ,n — TO — 1. 



tt 



(3.136) 



(3.14) 



3.3 Relation with the Hamiltonian reduction 



Theorem 3.1 provides an alternative definition of the space 6, and of the relative system of canonical 
coordinates: 

i = 1, . . . ,m 



<Ji = Qi = w*^* ^' 






5.1 



P^+. = 



dA 



Z-1, 
i= 1 



1=1, 



,n — m 



(3.15) 



We consider now the Hamiltonian (— Q) of the reduced f-flow, defined in Theorem 2.1 



SI 



= A-^K(g,), = A-^-— ^ 



Xi) 



^(5u(*+i)(t) 



At the end of the Chapter 2 we noticed how this expression looks very similar to a Legendre transform, but 
it is not; here we will show that actually the Legendre transform of the Lagrangian A gives the Hamiltonian 
Q, where Q is written in the coordinate system relative to A. 
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Firstly we rewrite Q in the coordinate system (3.15) 



-A + 


>; 




i=0 




m-1 


-A + 


>: 






,w 



i=0 
n — m — 1 n 



E 



61 \ au(^-i) 



^ '' j=m+l ^ ^ '' 



« \d^^^ (,) 



,w 






Using Lemma 3.1 we get 



m — 1 f. -f n—m~\ r^ » 

^=0 H 



E 



,w 

i=0 ^"tt 



= -A + ^p,( 



9ijt) 



for A(x, t, Qi, {qi)t, ..., (q„-„i)t)- 



(3.16) 



3.4 Concluding remarks 

The case considered in Theorem 3.1 is the more general one. Indeed, for {a — l)m < n < am, the Euler- 
Lagrange equation defining 6: 

^ • 



du{x) 
into the new "coordinates", is a differential equation in u, . . . , u^"^^^\ . . . , u("^™'"). 

The Lagrangian L transforms into 

L{x, t,u,u,,..., u^^^-^\uu . . . , u'^r'^\utu . . . , m("-™'"-i)) 
and we can define the new Lagrangian 

A/ J- (m-l) (m-l) (n-m-l.Q)\ 

A(a;,t, u,?/^:,...,^^ >,uu...,ul ',Utt,---,u^ '). 

The proof of the theorem is the same, one has only to consider the identities 

d f'dL\ _ d fdk'' 



dui^^l^) \dt J dui^'f^) V dx 
for z = 1, . . . , TO — 1 and /3 = 1, . . . , a. 
In particular, the manifold & is defined by 





(5M(™-i)(t) 
and it naturally carries the canonical system of coordinates 






u(*-i^^) i==l,...,TO; /3 = 0, ...,a-2 
^u'-'-^'"-^^ i^l,...,n- {a-l)m 

i — 1, . . . ,n 



(3.17) 



(3.18) 
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In the following wc will consider in detail the case a = 1, i.e. n < m, which occurs in the applications we 
are interested in (see next chapter). In this case L and L coincide and 



dL dL dL dL 



dt dt du * dui") 

The new Lagrangian is 



,(") 



"* + •••+ a^<- (3-19) 

A(a;, t, u,u^,..., m*""^', (u)t, ..., Uj""^') 



with 

dL dA dA dA dA ,. dA dA ^ ,^ ^„, 



H 



In this case Lemma 3.1 reduces to the following recurrence relation: 

SI dA 











Sui'1{x) 


du 


f 


and the proof is based on 


the 


identity 












d 


UL\ 


d 


(dA 








a«« 


Kdt)" 


dui'^ 


\ dx 


observing that, 


for i > 1, 








' ( 


dL\ 



i = l,...,n (3.21) 



— i^O,...,n-l, (3.22) 



dL 



du[''> \dt J du^^) ' 



In fact, L does not depend onu). , and 



d fdA\ d f dA\ dA 



du^'^KdxJ dx\Qu^)J auf-i)' 

In particular the first step, i ~ n, follow directly from the fact that the only dependence of Uj in both 
(3.19) and (3.20) is the one explicitly shown, so that 

^^-^^ (3 23) 

On the other hand, from (3.22) for the index i = 0, one obtains the fundamental relation 

SI _ dA du^''^'> 
Su{x) ~ du("^-^) dut ' 

o (m) 

and, since ^S — is always nonzero, the condition that defines the submanifold & is 

^^ 0. 



(3.24) 



5y(m-l) 

The relative system of canonical coordinates is given by 






du 

For « = 1, . . . , n. 
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The reduced t-flow is Lagrangian, with Lagrangian A. 
Indeed, from the identity 



du^) \dt J du(') \ dx 
and using the Lemma, one obtains on the subspace & 



d fdL\ d fdK\ 

-—] z = l,...,m~ 1, 



^^ ==0 i = n, . . . , m - 1 



which is the Euler-Lagrange equation for the Lagrangian A(a;, t,u,Ux, ■ ■ ■ ,u ,ut,---, u^ ) 



3.5 Example: KdV with tr fixed 

We wiU give below an example of how does Theorem 3.1 works for the first non trivial case, n — m. We 
study the Lagrangian reduction of the KdV equation 

ut = 6uux - Uxxx (3.25) 

on the stationary manifold of the ty-flow. 

The Lagrangian density of the ty-flow, reduced to the normal form (here I mean that L does not 
contains total derivatives), depends on the a;-derivatives oi u{x, t) up to order n = 3, and has the expression 

L = 7u^ + ibu^ul + 7umL + \{-^^^^f- (3-26) 

The submanifold S of the stationary points defined by the Euler-Lagrange equation for L gives the n + m = Q 
derivative in terms of the first five, explicitly 

u(*^' = 14wi('') + 2d,UxUxxx - 70u'^Uxx + 21mL - 70uul + 35m''. (3.27) 

From the relation 

dL _ dA 

dt dx 
one construct the Lagrangian A(u, u^, ■ ■ ■ , u^^'). By direct calculation 

A ^ - UV^ + i(.(^))^ - UuUxxU^^^ + lOuiui^yf + UUxUxxUxx.+ 

~70u'^UxU^^^ - {u^^^f + nu^{uxxf + TOuuluxx ~ 35u\xx+ 

35 

-—ul + 280u\l + 35u^. (3.28) 

The evolution equation (3.27) is the definition of Uxxx in terms of (u, Ux, Uxx, Ut), explicitly: 

Uxxx = Quux - Ut- 
Differentiating this relation in x one obtains 

u(4) ^ 6uUxx + Gul - Uxt 

u(5) — ISUxUxx + 36u^ — 6uUt — Uxxt 

u(6) =; 18u^^ + ISOUU^ + "i&U^Uxx - "iOUxUt - I2uuxt + Utt 
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which is a map from the "coordinates" 

u, Ux, Uxx, M*-^^ , u''^^ , u^^-* , u*-®-* , • ■ • 
into 

U,Ux,Uxx,Ut,Uxt,Uxxt,Utt, ■ ■ ■ 

The Lagrangian L depends on u,Ux,Uxx,Uxxx', in the new "coordinates" 

L{u, Ux,Uxx,ut) = 7u^ + bSu^ul + 7uul^ - QuUxUt + -u^ 
Its derivative 4t looks hke 

at 

dL dL dL dL dL dL 

-77 = ^TT + T^"t + ^ ^xt + ^ Uxxt + -^—Utt- 

dt ot du oux ouxx dut 

And there exist a functional A depending on u,Ux,Uxx,Ut,Uxt,Uxxt, explicitly 

A =35m^ + Au^ul + -ul - (Suluxt + -«'*+ 

-'i^U^Uxx - ^UU^^Uxx + SuUxtUxx + 11w^mL + 



-ul^ + GuUxUxxt + 22u^UxUt + AUxUxxUt - UxxtUt + 4UU( 



such that 



where 



dL 
dt 


= 


dA 
dx 


dA 
dx 


aA 

+ r. Ux 

du 


+ 


dA 

dux 


Uxx + 


dA 

dUxx 


^^(3) 




+ 


dA 

dut 
dL 

dt 


= 


dk 

^ Uxxt 
OUxt 


: + 


dk 

dUxxt 


,,(3) 







Ij A 9 9 

—— = SSu + lOGuu-. + Tu-.-. — QuxUt 
du 

— — = lOGu^Ux — 6uut 
dux 

dL 

= lAuUxx 

CfUxx 

dL 

- — = -Quux + ut 
dut 



and 



dA _ 
dx 



dA 
du 



210"^ + 12u u^ — lAQu Uxx — 2u^Uxx + ^UxtUxx + 22uu^^ + GuxUxxt + 4:4:UUxUt + 4mj 
dk o o 2 

:=: 8u Ux + 2u^ - l2UxUxt - AuUxUxx + GuUxxt + 22u Ut + ^UxxUt 

dux 
— = —35m — 2uu^ + %uuxt + 22u Uxx — Su^.^. + AuxUt 

OUxx 

dk 

—— = 22u'^Ux + 'iUxUxx - Uxxt + 8uut 

dut 

dk 2 

= -6m„ + Uxt + OUUxx 



duxt 
dk 



6uux — Ut 



dUx: 
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Lemma 3.2 states that the condition -?— = 0, which defines the submanifold 6, is equivalent to the condition 



dA d / dA 



duxx dt \ dux: 

explicitly: 

utt — 35u + 2'U'u^ — 22-^ Uxx + Suj-j, + 2uxUt — 2uuxt- 

The Euler-Lagrange equation for A reads 

— (- ) = Utt - 35u^ - 2uul + 22u'^Uxx - 3ul^ - 2uxUt + 2uuxt = 



duxx dt duxxt ' 

dA d , dA N -J ■=! ? 

7) If ^7) ' ^ ^"^^ + 2^2 - 4:UUxUxx - 2uUxxt + 22m Ut - '^UxxUt - Uxtt = 

dA d , dA N c 9 9 q 9 9 

'du ^ di^du' ^ "^ ^ 12m^< - UOu^xx ~ 2u>^^ + 4uxtUxx + 22uui^+ 

+ 2uxUxxt - 4Mf - 22u^Uxt + Uxxtt ~ Suutt = 0. (3.29) 



In this case L is nondegenerate, so that on 6 we can define the system of canonical coordinates 

fq, =«(*-!) i==l,2,3 
lK = ^7#T) *- 1,2,3 

which reads 

qi =u 

<72 = Ux 

<73 = Uxx 

Pl ^ 22u^Ux + 4:UxUxx - Uxxt + SwMt 

P2 = -6m^ + Uxt + ?>UUxx 

P3 = 6uux - Ut 



We will now solve the problem from the Hamiltonian point of view: starting from L and following Theorem 
1.1 one construct the canonical coordinates {pi,qi} on 6: 



qi 


= u 


12 


= Ux 


<?3 


^ Uxx 




51 


Vi 


SUx 




SI 


P-2 


6Uxx 




51 


Pa 


SUxxx 



70U Ux — 14:UxUxx — ^4:UUxxX + « 

= lAuUxx — M 



and the Hamiltonian function 

3 

1=1 
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By direct calculation one obtains 



35^" - UOu-'u 



3„,2 



35 



2 "= 



M„ — 35u Uxx + 70uu„Ux 



lu u^^ — u: 



XX XX 



OTrti ^CC 3?^^ I 



- ISUxUxxUxxx IQuul^x + (iulu^'^ + GuUxxU^'^ - ^l^^^))^ - 6uUxu('^ + UxxxU^'' 



and in canonical coordinates 
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Q - 35gf + 280g3g2 _ _q4 _ 35^4^^ ^ TOgig^gg - 21g2g2 + 

- ^3 - &qiq2Pi - Qqlp2 + 891(73^2 - -^pI - 70qlq2P3 - 4^293^3 + PlP3 + MipI- 
The corresponding Hamiltonian system reads 
qi = 6qiq2 - Ps 

92 = 6(72 -8^193 +P2 

93 = 709192 + 49293 - Pi - 891P3 

Pi = 21Qql + 8409?9| - 1409^93 + 709^93 - 42qiqj - 692P1 + 893P2 - I4O9192P3 + 4p^ ' 
P2 = 5609^92 - 7O92 + I4O919293 - 691P1 - 12g2P2 - 709^^3 - 493^3 
P3 = -35qt + 70qiql - 42qfq'^ - 3ql + 8qiP2 - M2P3 

Rewriting this system in coordinates {pi,qi} one obtains exactly (3.29). 



3. A Appendix 



Proof of Lemma 3.1: wc prove the Lemma in two parts: 
• firstly we prove the relation 



E 






^\SuUHx)J Ouf^ 



i — 0, . . . , n — ?Ti — 1. 



(a.l) 



For convenience we can explicitly rewrite eq. (3.8a), using (3.9): 

Bt J^ ^ UuO) J dt 

j—7n 






i=0 
n — rn 



n 

E E 






,w 



i=0 '-j=m+i ^ ' '^"t 



,W 



(a.2) 
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Notice that the arguments of in the square brackets depend on u and its a:;-derivatives upon the order m — 1 
and on ut and its x-derivatives upon the order n — m, so that the dependence on Uj for n — m + l < i < ra—1 
and on Ujj , for every j is only the exphcit one. Analogously 



dx 



dk dk au(™) 
a^ ^ au("'-i) dt 

dk au(™) ■ 



E 



dA 



du(i-^) 



,w^ 



5u(™-i) du 

dk du^"^'> 
5^™-!) dut 



Ut 



-utt 



E 

n— m 



dk du^""^ 



dk 



duf^ du^r~^^ du(^ 



,{i) 



dA 



M) 



The i-th step of (a.l) is obtained from the obvious identity 



d fdL 



duif\dtj du^}\dx 



d fdA\ 

-r\ « = 1,. 



Indeed, from (a. 2) and (a. 3), it follows that 

d f'dL 



(i) V dt 



and 



du. 



d fdA 

dx 



E 



dL \ du(J^ 






dA 



dA 



In particular, at the first step, i ^ n — m one obtains the basic relation 



dA 



dL \ du^"^ 



61 



du'-^' 



Substituting (a. 6) into the further step of the recurrence, one finds 

dA / 'sT \ 9u("-2) / JT \ du^"-^^ 



and so on. This gives relation (a. 2). 



(a.3) 



(a.4) 



(a.5) 



(a.6) 



• The second step is the proof of the relation 

dk / IT \ A/«^\ du^^-''> 



E 



UUf. \ \ 1/ 2=ra+\ ^ ^ >' 



i ^ n — m, . . . ,m — 1, 



(a.7) 



which is a part of (a.l), indeed, for i > n — m, the partial derivative ^^i) vanishes. Very much as in the 
previous case, equation (a.7) follows from the identity 



dui: 



d fdL\ _ d fdA\ 



dutt \ dt J dutt \dx J 
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Using (a. 2) one can rewrite the left hand side as 

d flh 



dutt V dt 



E 

j=m 



dL \ 9u(j) 



duii) du 



On the other hand, from (a. 3), one has 

d fdk 



d f dK 



dk 9u('") 
dutt\dx ) dx\duu) ' du]"^^^^ dut 



where 



d f dk 



E (; 



IT , 9u(J-i) 



dxKdutt ) dx\ ^—' ^Su'^^Hx)' dut 
We develop the right hand side, recalling (3.6b): 



dx \ dut 



dut 



obtaining 



d f dk 



dx \ dut 



E 

j—7n-]-l 

d 



61 \ 9u(j'-i) 



5y(m-l) Q^^ 



SI \ du^^'> 



61 



Q^(m) Q^{m) 



dx\6u^J'>{x) J dut \6u^i'>{x)J dut \6u'-i'> (x) J du^'^^^'l dut 



61 



_dx\6u'y'^'+'^){x) 

'JT 



9u(™) 



E 



E 

d 



dut 



/ 61 



\ 6u(i) (x) J dx\ (5^0+1) (a;) 



61 



(5u(j)(a;)y 5m("'-i) dut 



61 



dx\6u("'+'^^x) 



dut 



E 



j—m+l 

Inserting in (a. 9) and equating it to (a. 8) we obtain 



dL \ 9u(j) 



.5/ 



du(^'>J dut \6u(^\x)J du("'-^'> dut 



du^^l 
dut 



Q^(m) g^im) 



f dL 



_Vau(™)/ dx \6u('^+^'> (x) 



51 



61 



Q^(m) Q^(m) 



du^""^ _ dk du^""^ 
dut ~ au(™-i) dut \6u<-i^x)J au^"-!) dut 



But the term ".. is nonzero by definition, so that 



dut 



dk 



61 



E 



61 \ du^i-^'^ 



CUj \ \ )/ j=m+l ^ ^ ' '' 

which is the first step of the recurrence (a. 7), and so on. 



(a.8) 



(a.9) 



• Finally, since (a.l) for i > n — m coincides with (a. 7), it remains to prove it for i < n — m. These 
relations can be obtained from (a. 2) and (a. 7) together with the identity 



d fdL 



dui'^\dtj dui^'^Kd^ 



d fdk\ 

-;- , » = 1, 
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Indeed, starting from the index i = n — m and using (3.9), one may write 
d f'dL\ d f dl \ dl 



d f dt au(") \ f 'di 






On the other hand 



d fdk 



dui"-"'^\dxj dx\Qj^---^ 



dA 



dA 



Performing the same steps as in the previous case, one obtains 



(n— m— 1) 



+ 



dA 



du 



(m-l)- 



dA 



■(- 



dA 



61 



n , 



(5/ 



which is the first recursive step of (a.l). 



(a.lO) 



Q.E.D. 



Proof of Lemma 3.2: We prove the Lemma by mean of the equivalence 



61 

Su{x) 



dA 



gy(m-l) rftV(9u'™~^^ 



dA 



dut 



d fdL\ 



d fdA\ 



which foUows from the identity 



dut \dt J dut \dx J 
Indeed, expanding, one has 

d fdL\ _ d dl dl _ 
dut \ dt J dt dut du 

~Jt^ \dui3)j dut 



dL 
du 



E 



dL \ du(^^^ 



d 
dt 



dA 



dA du^""^ 



dx\duttj du^r ^^ ^"* 



du^j") ) du 
'd^J ^ ^ \du(^)) du 



where the last equality follows from the equivalence of (a. 8) and (a. 9). Expanding the right hand side 



d fdL 
dut \ dt 



On the other hand 
d /dA 
dut V dx 



d 
dx 



dA 
dt \dutt 



dA 



du 



d I dA 
dx \dut 

d_fdP^ 
dx \dut 



dL\ ^ 



]=m 



dtKdui""-'^ 



dL \ du'^'-i^ 



du^""'' 
dut 



dA d fdu("''> 



du'r'^'>dt\ dut 



du'^J'i ) du 



(a.ll) 



dA du^""^ 
gy(m-i) g^^ 

dA du^"''> 
9u(™^i) dut 



dA du^-^^ 



dA d /du^""^ 



ay(™-i) dt \ dut 



E 



SI 



auO-1) 



Su("''>(x)J ^—' \3u^i){x)) du^"^-^) 
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dA d /du'^"'^ 



Qui^) 



du au*""^' dt \ dut 



(a.l2) 



Comparing (a. 11) and (a. 12) one obtains 



d 
dx 



dA\ d / dA 



duf J dt\ dutt 
du I ' ^-^ V 9u(j) / du 



dA 



dA 



-(l^) + fl 



'IT 



But Lemma 3.1 states that 



dA\ d f dA 



so that 



d_ 
dx 



dut J dt \dutt 
dA\ d f dA 



dut I dt \dutt 



(5M(™)(a;) 



SI 

Sux{x) 



SI 



dut 

=m+l ^ ^ ' ' 



j—'m-\-l 



E 



du 



J— m+l 



SI \ du^^-^^ 



^\Su(i'>{x)J du 



dx \ Sux (x) J dx 

n 

E 



E 



SI 



V<5m(J')(x) 
f 'TT \ d fdu^^-^^ 



j=m+ 

d / 'sT \ d 



\ Su'^i^ (x) J dx \ du 

n 

E 



dx \ Sux (x) J dx 



SI 



dvi^ 
du 



du 



j—m-]-! 






^ \Su(^Hx)J\ du au("-i) du 



Substituting in (a. 13) we get 



dA 



Qu{m-1) J dt\Qyj^"^-^) 



dA 



dut 



E 

j=m 



dL\ d r SI 

du J dx \Sux{x) 

aL \ d 



duii) J dx\Su''J+^){x) 



SI 



du'^'-i^ 
du 



J— m+l 



SI 



du'^^^ du^^-^^ au("-i) 



^ \5u(3)(x))\ du dui"^-'^) du 



(a.l3) 



SI 



(5m(™)(x) 



E 



SI \ du^i-^^i 



^ \5u(i){x)J du^-^-^'^ 



du 



All the terms cancels but 



dA 



du^rn-l) J dt\Qyl''^-'^) 



dA 



du^'^ 



dut 



SI 



dL 



du J dx \ Sux {x) 



Q.E.D. 



Proof of Lemma 3.3: Relation (3.13) follows from the identities 

i == 1 , . . . , TO — 1 , 



d fdL\ _ d fdA'^ 



du'y') \dt J dui^ \dx J 
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(a.l4) 



and 



d fdL\ _ d (dA\^ 



du^^\dtj dui'^\dxj 



2 = 1 , . . . , JTl — 1 . 



Starting from (a. 14), wc can write 



d fdL\ _dfdL\_d^f 'dL \ 9m(j) d f 'dL 



E 



9u(*) \dt ) dt Vau(») ) dt ^ \duio) J dui') dt \du(') J ' 



]=m 



and 



which give 



d /dA\ d f dA 



dA 



dk au(™) dk d fdu<-"'^ 



du^'^ydx) dx\du(^^ ) du^'--'^') i9u(™-i) Su^*) Qu^"^-^^ dt \ du^'^ J ' 



d_ 
di 



E 



dL \ 9u(j) / dL 



duO) J du('') \ dui^) 



d f dA 



dA 



dk au(™) dk d fdu^""^ 



dx\du(^ J dui^-^) au(™-i) auW Qyi^'^-^)dt\du(^ 



(a.l5) 



(a.l6) 



On the other hand, in (a. 15) 

d fdL\ dL d f dL 



dufydtjdu'^^ dt\duf) 



dL 



and 



which give 



d fdA\ d f dA\ dA 



E 



f dL \ du^^^ d f dt 



Kdu'^j) du^"^ dt\oJ'^ 



dA du^-"'') 



dL 



E 



f 'dL \ du^^^ d f dL \ d f dA \ , dA dA du^"''^ 



Performing the derivative w.r.t. t and substituting into (a. 16) one gets 





- 


d^ f dt 

dt' U^,W 


\ d d f 
J dx dt\ 


dA\ d dA 
dufU ' dtd^f-^ 


) + 


d 
dt 


■ dA 

Idul 


du^'^y 

1) dui^ 


= 




_ d / dk 
dx\du(') 


\ dk 


dk du^"''> 


dk d 
^ du'T''^ dt 


^au(")\ 




ay(m-l) gy^(z) 


, dui^ J ' 


which gives 










d^ f dL\ 
dtAdui'^J 


d 
dx 


'( dk\ d f dA Y 


+ 


" dk 

dui'-^) 


d dA ' 
dtduf^l 




+ 


" dk 


d aA 








5y(m-l) 





The left hand side of (a. 17) is zero ii i > n ~ m 
li i < n — m, 

dL 



du 



(«) 



, ^ \d^)^^ 

k=m+i ^ ^ '-'"■t 



dA 



dA 



9u^r'^ '^^ ^duif 



(a.l7) 



(a.l8) 
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Indeed, using (a.l), and relation (3.6b), which we rewrite here below. 



dx\duf^J dui'^ 5w(»-i) 



one obtains 














= 


n 

- >: 


' d 1 SI Y 




\ aM(*+™-i) 


dxlv<5MW(a;)j 


)y 9z.f-^) 


Z^ l^y(fe) )^ w + 


" d / 61 Y 






_rfxUw(™+*+iHa;)j. 













where the last identity follows from the fact that 










li k ~ I < m. In (3.6), this implies 



Qy^(k+m) g^{k+m+j) 



du. 



(fc) 



du'^''+'^ 



Finally, substituting (a. 18) into (a. 17) gives (3.13) 



dk 



dA 



9u(-i) dt\dut'^J dt^\gu[l-^) 



dA 



dA 



5y(m-i) (H Vau*""^' 



dA 



au(™) d 



dk d f dk 



du(^) dt\Quf) dt^Kdv!^ 



4. Applications to Painleve equations 



dA 



Q.E.D. 



In this Section we study some applications of Theorem 2.1 and Theorem 3.1. to show how the finite 
dimensional Hamiltonian structure of Painleve equations comes from an infinite dimensional structure via 
the above procedure. 

4.1 PI as scaling reduction of KdV 

At the beginning we study the problem following the Hamiltonian scheme, then we will apply the 
framework of Theorem 3.1. 



We consider the KdV equation 



Ut ~ 6uUx — Uxxx- 
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(4.1) 



i.e. the t — i^-flow in the KdV hierarchy (1.2); it admits the nondegenerate scaling symmetry 

f u^ 

I^^^ := [u^ + -^ + 2ux + Uu^)dx, (4.2) 

which depends on x,u,Ux,t. We note that L = [L^ + AxL-^ + 12^^^], where I_^ = j L_^dx — J ^^^^dx and 
I„ = J Lgdx = J —j^dx are the first Hamihonians of the KdV hierarchy. 

Theorem 2.1 states that the t-flow is Hamihonian on the manifold 6 of the stationary points of the symmetry, 
i.e. 6 is the 2-dimensional manifold of the solutions of the Euler-Lagrange equation 

A J 

= Ua;x - 3^2 - 2x - 12tu = 0. (4.3) 



6u(x 

It is invariant under the t th flow and it naturally carries the system of canonical coordinates 

q = u 

P = S — = Ut 

Notice that the identities 

hold, where H is the generalized Lcgendre transform of L: 



(4.3) 



51 

H = -L + - — Ux. 

OUx 

The first of identities (4.4) allows us to express the higher derivatives u*^™^ for to > 2 in terms of x,t, p, q 
and pw with I — I, ... ,171 ~ 2 + I. 

On & px + 4^ = 0, and the system (4.4) reduces to the canonical Hamiltonian system 

Px ^ 'iq^ + 2x + 12tq 
Qx =P, 

for the Hamiltonian function 

H = -L + ul ^ ^ - q^ - 2qx ~ 6tq^ (4.5) 

giving the reduced a:;-flow. This system is equivalent to the second order ODE in the variable q : 

q = -2,q^ -2x- I2tq. (4.6) 

The space 6 is the set of the stationary points of the scaling symmetry (4.2); this means that & carries a 
"natural" system of canonical coordinates {w^, tt^}, given by the self-similar function of u, i.e; combinations 
of u,x,t in the variable z{x,t) invariant w.r.t. the scaling. We will call them scaling coordinates. In this 
case 

TT = 2p 

with z ^ X — 6t^. 

In terms of the scaling coordinates the system reads 

{dn _ af) 
dz dw 

dw ^ dS) 

dz dTT 
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for the Hamiltonian 

iO = — - SW^ - Awz + Si'' + Atz. 
8 

The system is equivalent to the ODE: 

w" ^ 6u;^ + z, (4.7) 

that is exactly Painleve I. The Hamiltonian fy differs from the usual PI Hamiltonian for the terms in z, t 
that do not enter in the Hamiltonian system. 

We now construct the time dependent Hamiltonian function (— Q), that is the reduction on 6 of 

where p, q are expressed in terms of u, u^, and A(a;, t, u, u^, Uxxi Uxxx), calculated from 

dL _ dK 
dt dx 

has the form 

9 1 

A = 6i(4u^ - 2uUxx + ul) + 2x{3u'^ - Uxx) + -ZU^ + -^ulx + '^'^x - Su'^Uxx + Quul - UxUxxx- 

By direct calculation one obtains 

y2 y2 g 

Q = 12t{2u^ + ^ ~ ""^^) + ^ " Su^uxx + 2" + 2ux + 2x{3u^ - Uxx), (4.8) 

This reduces on 6 to 

„2 

Q = 12t(y -q^ - 6tq^ - 2xq) + 2p - 2x^ (4.9) 

Theorem 2.1 states that (— Q) is the Hamiltonian for the reduced t-flow,i.e., in terms of p and q 

<^ ^P - (4 10) 

[ p = -m {3q^ + 2x + 12tq) = ^, 

Notice that system (4.10), written in terms of the scaling coordinates w and z, gives the same Painleve I. 

Remark: In this case the evolution equation is Hamiltonian and it can be written in the form 

u,^{uixlh} = - — , 



where I^ = J L^dx with density 






On the other hand the scaling symmetry defines the stationary flow 

du 

-— = 12tux + ut + 2 = 0. 

ds 

which is Hamiltonian: 

— = {u{x),I} = — — ==0. 
ds dx ou 
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The s-flow and the t-flow commute, but the Hamiltonian generating the scahng depends exphcitly on 
the time t, so that the relation 

dl, , dl, , 

—^ ^ \i I \4 — <;i = 
dt ^'(=)''i^+ dt 

holds. Hence we have an alternative way to define the reduced Hamiltonian Q, following [BN]: 

d dL dl d dh 

dx dt du dx du 

In this case the relation (4.2) follows as a consequence. 



System (4.10) i.e. the reduction of the i-flow on (5, can be obtained from the Lagrangian point of view; 
indeed one can consider the evolution equation (4.1) as the definition of tixxx in terms of {u^UxTUxxT'U't)^ 
explicitly: 

Uxxx = QuUx — lit- 

Differentiating this relation in x one obtains 



l'-^^ = QuUxx + dul - Uxt 

i(5) — ISUxUxx + 36u^ — 6uUt — Uxxt 



which is a map from the "coordinates" 

U,Ux,Uxx,U^^\u^^\u^''\u^''\... 

into 

U,Ux,Uxx,Ut, Uxt, Uxxt, Utt, ■ ■ ■ 

The Lagrangian L depends on u, u^ and hence its derivative -^ looks like 

dL _ dL dL dL 

dt dt du dux 

And there exist a functional A depending on u, Ux, Uxx, ut such that 

dL dk dk dk dk dk ,3^ dk 

dt dx dx du dUx duxx dut 

Here, in terms of the new coordinates 

9 4 1 
A = 6t(4u^ - 2uuxx + ul) + 2x(3u^ - Uxx) + ^w + ^"L + "^"^x - 'iu^Uxx + UxUt, 

L = L, and 



dL 

dt 


= 6^2 


dL 

dv 


= iu^ + 2x + 12tu 


dL 

dux 


= Ux 
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9A ^ 2 „ 

ox 

dk , „ 

—- = 72tu^ - I2tuxx + 12a;M + IM^ - Quu^ 

ou 

5A _ 

-— = \2tux + ut + 2 
dUx 

dk 

— —12tu + Uxx — 2x — 3m 



dk 



dut 
The condition 4^ = 0, that defines the subnianifold S, is equivalent to the condition 

'^ 



O^xx 

Hence we have an alternative definition of the space 6, and an alternative way to defines the canonical 

coordinates: 

q ~ u 

Theorem 3.1 states that the reduced t-flow is Lagrangian ,with Lagrangian A, in this case it is easy to verify 
it, indeed, on 6, 



-12tu + Ujjj. — 2a; — "iu^ = 



dk 
^ = I2tux + ut + 2 = Q 

dk d / dA 



. ^ - fti§^) = 72tu2 - UtUxx + I2xu + 18^3 - QuUxx - Uxt = 0, 

where the first equation is the definition of the submanifold 6 itself, the other two reproduces (4.9), indeed 
they can be rewritten as 

ut = —12tux — 2 

u^t = -12t(3u2 + 2.x + I2tu) 



4.2 PII as scaling reduction of mKdV 

One can repeat the same procedure as in section 4.1 starting from the mKdV equation 

Ut — 6u Ux — u^^^. (4.11) 

It admits the nondegenerate scaling symmetry 

(^t{u^ + ul) + ^)dx, (4.12) 

2 

which depends on x, u, Ux, t. We notice that L — 3tL^ -\- ^^. 

Here © is the 2-dimensional manifold of the solutions of the Eulcr-Lagrange equation 

Uxx - iri^^^^ + "^) =" 0. (4.13) 



Su{x) 'M 
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It naturally carries the system of canonical coordinates: 



q ~ u 

SUj: 



P^ Ir ^ 3tMj 



As in the previous case we read the Euler-Lagrange equation as a reduced x~flow with Hamiltonian 

H = '---tq^--q^. (4.14) 

where H is the generalized Legendre transform of L: 

SI , 

H = -L+ - — Ux = -L + Stui.. 
Sux 

The system is equivalent to the second order ODE in the variable q : 

. 1 
Qxx = 2g + —qx. 

The scaling coordinates are now 

w = (3i)*<7 

(3t)3 

in the variable z — — ^, and the system transforms into 

(3t)3 

{dn _ d^ 
dz dw 

dw _ df) 

dz d-TT 

for the Hamiltonian 

2^ '2 

The system is equivalent to the ODE: 

w" = 2w^ + zw. (4.15) 

that is exactly Painleve II. 
We now construct the time dependent Hamiltonian function (— Q), that is the reduction on & of 

-Q^-A + 4, 

at 
where 

A = t{6u^ - 6u^Uxx + ISu'^ul + -ul^ - SUxUxxx) + x{-U^ + -ul - UUxx) + UUx. 

By direct calculation one obtains 

13 1 

Q ^ 6t{u^ - u^u^x + -uU + x{-u^ - uu^^ + -ul) + uu^, (4.16) 



which on & reduces to 



f<l^-'''+s'"-r'^-^ '"" 



and is the Hamiltonian for the reduced i-flow . In fact 



y 3t W dp 

• _ p_ _ ££f^ _ 2 3 _ dQ 
P ~ 3t 3t ^^ -^ ~ dq' 
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(4.18) 



Notice that also the system (4.18), written in w and z, gives Painleve II. 

Remark: The evolution equation is Haniiltonian and can be written in the form 

d SI 

where /^ = J L-^dx with density 

On the other hand the scaling symmetry defines the Hamiltonian stationary flow 

du ^ d SI 

— — — XUx + 6tut + u — r — = 0. 

ds dx Qu 



We now deduce system (4.18) from the Lagrangian point of view, reading the evolution equation (4.11) as 
the definition of Uxxx in terms of (m, u^, u^x, Ut), explicitly: 

Uxxx = Gu^Ux ~ Ut. 



Differentiating this relation in x one obtains 



U^^' = 6u^Uxx + 12uUx — Uxt 

U^^' — iQuUxUxx + iQu^Ux + 12ui^ — Qu^Ut — Uxxt 



which is a map from the "coordinates" 



into 



Here 



and 



U,Ux,Uxx,U^''\u^'^\u^''\u^^\ 



U,Ux,Uxx,Ut,Uxt, Uxxt ,Utt, 



A = t{6u^ - 6u^Uxx + -^ulx " iUxUxxx + iUxUt) + x{-u'^ + -U^. - UUxx) + UUx 

oL 3,4 2 \ 

— — = Dtw + UX 

au 
- — = Mux 

OUx 

^^ 3,4, 2^ 

T^ = -(U +Uj -UUxx 

dx 2 

dA t; 2 s 

— — = '36tu — 18tu Uxx + 6xu — xUxx + Ux 

ou 

dA 

— = XUx + U + itUt 

OUx 

dA ^ ,0 

AtUxT — btu — UX 



dux. 



dA 

—— = StUx 

dut 
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The condition t— = 0, that defines the subnianifold 6, is equivalent to the condition 

dA 



dz 



= StUxx — Qtu — ux = 



Hence we have an alternative definition of the space 6, and an alternative way to defines the canonical 

coordinates: 

q — u 

On 6 



dA 

9Mxx 

dA 
dux 



'^tuxx ~ 6tu^ — ux = 
■§^ = XUx + itut + u = Q 



where the first defines the subnianifold, the second one gives the motion of u and the third the motion of 
Ux, hence one can rewrite them as 

r itut = —XUx — u 

\ StUxt ^Ux- 2u^x - 3^, 

which coincides with (4.18). 

4.3 PHI as scaling reduction of Sine-Gordon 

A particular case of Painlevc III equation can be obtained as reduction of the Sine-Gordon equation 

(4.19) 
{Vt ^ Uxx - sm u ^ --^, 

via the scaling 



51^ 

vt = Uxx - sm u ^ -7:7, 



S -^ OV 

^ = x{uxx - sin u) + tvx + Ux = — j^ = 0, 



^^XV + tUx = ^=Q 

SI, 



(4.20) 



where I^ — J L^dx, I — J Ldx, with the Hamiltonians 

1 



and 



w.r.t. the Poisson bracket 



L, = -{v +uJ-cos u 



L = —{v + Uj.) — XCOS U + tvUx — xL^ + tvUx 



iFr\^ f( ^^ ^^ ^^ ^^ ^ d 

^'^ J^5u{x)5v{x) 5v{x)5u{xy 
The scaling reduction equation means 



'^^(e, _ '5^<e, _ 



5u{x) 5v{x) 
which defines the subnianifold S: 

XV 

Ux - Y' 

with the canonical coordinates 

p = XUx +tv — ^^ V 

q = u 
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(4.21) 



in 6. 

The equation defining 6 can be written as an Hamiltonian system in canonical form describing tlie 
reduced x-flow: 

/ W^ = -f 
I i^)^ = f 

where 

2 

^=1^2^+^ cos g. (4.22) 

In terms of the scaling coordinates 

in the variable z — ^ „ , the Hamiltonian system transform into 

{dn _ of) 
dz dw 

dw __ df) 

dz Ott 

for the Hamiltonian 

l7r2 
Sj ~ cos w. 

4 z 

The system is equivalent to Painleve IH: 

2zw + 2w — sin w = 0. 



Let us now construct the time dependent Hamiltonian function (— Q), that is the reduction on 6 of 
where 

A = XUxV + t( — {v +u^)—cosuj. 



By direct calculation one obtains 



which on & reduces to 



t{^iv^ -ul)-cosu) 



-'^l^^-'"'"'^- 



This is the Hamiltonian for the reduced t-flow . In fact 



1„ - -^^P-- - _M 
p = V + xVx + tvt — t sm q — -rf- . 



LJ — (y — ~ LL^r — " t/ , o o — ~ -^ 

■i X -^ t'^—x^ op 



dq- 

Note that also the system (4.23), written in w and z, gives Painleve III. 

Remark : We now deduce system (4.23) from the Lagrangian point of view, reading the evolution equation 

(4.19) as the definition of v in terms of ut, explicitly: 

(V^Ut 

\ Uxx ^ vt+ sin u = utt + sin u 
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Differentiating this relation in x one obtains 

Vxx = Vtt — V COS U — Uttt — Ut COS u 
Uxxx = Uttt + UxCOSU 

which is a map from the "coordinates" 

lA, V ^ Ux^ Vx^ Uxx: ^xx^ • • • 



into 



Here 



and 



which give 



and 



U,Ux,Ut,Uxt,Utt,Uxtt,-- ■ 



A = xUxUt + t(^ — {ut + u^) — cosu) 



L — tuxUt +x(-(uf + u^) — cosu) 



dL 1,2 2 N 

dl 

— — — — sm u 

ou 

dl 

= tUt + XUx 

OUx 

dl 

= tUx + XUt 

out 



9^ 1/2 2 ^ 

-Q^ = 2^ +u^)-cosu 

dk 

— — = — tsmu 

ou 

dk 

= XUt + tUx 

OUx 

dk 

= XUx + tUt 

dut 



The condition 4^ = 0, that defines the submanifold S, is equivalent to the condition 

dk 

= XUt +tUx=0 

OUx 

Hence we have an alternative definition of the space 6, and an alternative way to defines the canonical 

coordinates: 

q = u 
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On 6 

I li = a^"* + ^"- = 

where the first defines the submanifold, and the second 

—t sin u — xu^t -~ Ut — tutt — 
coincides with (4.23). 

5. Self— similar solutions of nwaves equation and Hamiltonian MPDEs 

5.1 n— waves equations and their symmetries 

Let us consider the equation 

ut-v^ - [u,v] = 0, (5.1) 

where 

u=[j,a] v—[j,b] a = diag {a^,....a"-) b ~ diag {b^, ....b"-) (5.2) 

and 7 is a function of x, t. 

Following [DS] it is possible to rewrite (5.1) as an infinite dimensional Hamiltonian system on the 
space 971 of functions of x with values in Mat{n, C) with the inner product 



(w,w) — Tt {u{x)v{x))dx. 
On the space ^ of functionals 

F= f{x,u,Ux,....u ) dx 
one can define V^F G 9Jl by 



d 
'de 
and the Poisson structure P with the Poisson bracket 



F{u + ew)\^^^ = (V„F,w) 



{i^,G}(^i) = (V„F,[V„G,^+«]) (5.3) 

The n-waves equation (5.1) is a Hamiltonian system w.r.t. this Poisson structure: 

ut = Pdh = [VJi, — +u] = [-V, ^ + "], (5.4) 

where 

Ii ~ I Lidx = / Tr {uv)dx 

that in components of 7 gives 

h= / X!X!( — ~ — u,kUki)dx^ ^^[{h-bk)ia^-akh,^-f,^,]dx- (5.5) 

I k I k 
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For (n = 3, M ~ —u) one can reduce to a particular case of P VI equation (see [D], where the 
Haniiltonian structure for this particular case of P VI is derived from the Hamiltonian structure of the 
n-waves equation.), imposing the scaling 

du „ ,. „, 

— — = tUt + XUx + M = 0. (&-6) 

as 

It admits the Hamiltonian form 

J.[V^,,,.4+„|.0 (5.7) 

where 

/(s) = / Ldx ^ -- Tr {tuv + xu^)dx = / X! X!^^ — 'Z 7;)'^ikUkidx 

J ^ J J ■ , ^k 0,1 Z 

I k 

and 

V^/(s) = —tv — xu, (5-8) 

or, in terms of 7: 

^(«) ^ '^'^ii^^^^k){ai- ak)t+ (a^- ak)'^x]j^^-f^^dx (5.9) 

i k 

We emphasize the fact that the t-flow and the s~flow commute, so that 

dL ^ 

By substituting: 

/'[Tr(V„/iK/(,),— +«]) + lTr{uv)]dx - 0. 

Then there exists a function Q {x, t, u, v) such that 

d uv d 

Tr{-v[-tv-xu,-+u] + -)^--Q^^^. 

By direct calculation (see Appendix 5. A) we obtain 

Qw = \Tr{xuv + tv^) - i ^[(a^ - a J (6^ - b^)x + (6^ - 6j't]7.,7,. 

As in the previous examples, Q,^ is the Hamiltonian for the reduced i-flow. We now describe this 
flow. 

We start by rewriting the system 

ut-v^- [u, u] = .g -^Qx 

tUt + XUx + U — 



in terms of 7, i.e. we solve 
under the condition 

This gives 



[it, a] = [7^, 5] + [[7, a], [7, 6]; 



t 1 

Ix = It 7- 

X X 



[yt,ax + tb] + [7, b] = [[7, ax], [7, b]] 
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but, because of the commutativity of b with itself, 

d 



[7, ax + tb] = [[7, ax + bt], [7, b]]. (5.11) 

at 



Then we identify 6 s with the space of matrices 

q ~ [7, ax + bt] = ux + vt, 
so that u = ad, , ^,,ad^^q, v — ad, , ^,,ad^^q, or, in terms of the matrix elements: 

(ax+bt) a ^' (ax+bt) b ^' ' 

9,, = [(a, - aj.x + (6^ - 6Ji]7,, 
On & the equation (5.11) has the Lax form 

qt = [q,v] = [g,ad(^^^^^)ad^iq] 
with the Hamiltonian function 

2 '^ ' 2 

This coincides with Q . 



1^, . 1^, 



^(t) = oTr(gw) = -Tr(a;uw + tw^). (5.12) 



One may change the role of x and t. This means that one considers the system (5.1) on the space of functions 
vit): 

Vx = \^vla,-i: +v], 

where Iq = J H^dt = — i J Ti{uv) dt and one integrates in the variable t. The scaling (5.8) can be read as 
an Hamiltonian equation 

dv , d ^ , , 

^^[V./(.),- + .]^0, (5.13) 

where 

/(,) = Ldt = -- / Tr {xuv + tv"^) dt (5.14) 

and 

V^,/(s) = —tv — xu (5.15) 

Commutativity of the flows is equivalent to 

dl, , 

in our case 

[Tr(V„/o[V„/(,), ^ + «]) + \Tr(uv)]dt - 0. 

Then there exists a function Q {x, t, u, v) such that 

, . d , uv^ d 

T,(_u[-tv-xu,-+v] + -)^^-Q^^^. 

By direct calculation (see Appendix 5. A) we obtain 

Q(,) =- ^Tr{tuv + xu^), (5.16) 
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(5.17) 



in components: 

Now we study the a:;-flow on the reduced manifold defined by the scahng equation 
the system (5.10) gives 

ut -v^ - [u, v] =0 

tVt + XVx + V ~ 0. 

In terms of 7 this becomes 

— [7, ax + th] = [[7, ax + bt], [7, a]], 
ax 

that is a Lax equation on S^: 

qx = [q, u] = [q, ad(_^^^^^jad;;"^q] (5.18) 

with Hamiltonian function 

•^(x) = 9Tr((7u) = -Tr(xu2 + tuv). (5.19) 

This coincides with Q . 

In fact one can rewrite the scaling as a zero-curvature equation in two ways: 

-^^q^ + [u,q]=0 (5.20) 

as 

and 

du , , , , 

— = q, + [v,q]=Q. (5.21) 

as 

Therefore one may rewrite them in terms of q as 

qx — \q, ad'' ad q] 

and 

qt — \q, ad'' ad.q]. 

5.2 Commuting time— dependent Hamiltonian flows on so(n) 
We can do exactly the same using the coordinates 

and the corresponding derivatives ^, with 

d y--^ d 

and 

d Y^ , "^ 

i ^ 

The starting equation is now 



d d 

-—Ui~—-Uk-[u^,Uk\^0 (5.22) 

at]^ at . 
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where 

u, = [7, Ei\ {ui)ki = l^l5^k - IkiSa 

and {Ei)ki — Sik^ki- We impose the scahng 

d ■r-^ d 

-f'^k = 2_^U-7-Uk + Uk=0 (5.23) 

i * 

For every k one can define, on the space ^^ of functionals 

/du d"^u. , 

■'^ ' ' dt' ' dt]^' 

with 

—F{uk + ew)\^^„ ^{V^^F,w), 

a Poisson structure P with the Poisson bracket 

{F, G}{uk) = (V„^ F, [V„^ G, ^ + «,]) 

The n-waves equation (5.22) is Hamihonian w.r.t. the Poisson structure P in ^k'- 

where 

Lidt = -- / Tr (uiUfe) difc- (5.25) 

On 5^fe we can reduce to P VI equation imposing the scaling (5.23), which admits the Hamiltonian 
form 

j 



I{s) = / Ldtk = ^^ / T^'f'^tjUjUkdtk = ^^ / Tr(^tjUjWfc + tkul)dtk (5.27) 

The commutativity of the flows is equivalent to 



in our case 

[Tr(V„^/4V„^/(s), — +Ufc]) + -Tr(uiUk)]dtk == - / dkQ^^^dtk- 

By direct calculation (following the scheme in Appendix 5. A) we obtain 

Q<M ^ ^Tr^t.u.^i, - J2^U - tj)7.,%. (5.28) 

j 3 

The scaling equation defines the submanifold Sg. One can consider on &a the system of coordinates 
given by the matrix elements of q: 

9= [l.^-tjEj] ^ [7,t/], 
i 
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where U is the diagonal matrix diag(fi, . . . , t„); exphcitly 

Qv = (t, - t^h.r (5-29) 

As in the previous cases, Q,. is the Hamihonian for the i^-flow on the reduced manifold. 

Starting now from ^^, i ^ k, we can reduce on the same submanifold & and construct the Hamiltonian 
function Q. . 

Indeed, the scaling (5.23) for every k produces on &s the Lax equation 

Qk = [q,Uk], (5.30) 

with Hamiltonian functions 

Hk = iMqu,) = lTrJ2t,u,u, = i^M^. (5.31) 

These coincide with the Q.^., constructed above. Observing that 7 = ad^J^q one can rewrite 

Uk = adj,^ad~^q. 

In the case q = — g eqs. (5.30) are the Monodromy Preserving Deformation equations for the linear 
differential operator 

that give Painleve VI, for n = 3, and the higher-order analogues, for n > 3. 

Remark: The first integrals of the MPDE (5.30) are given by the monodromy data of the operator A. The 
Poisson bracket on the space of the monodromy data has been computed in [Ug] . 



5. A Appendix 

Here we present the explicit calculations giving rise to equation (5.16): let us consider the following 
explicit expressions: 

/f = -- / Tr {uv)dx 
Vlt - -V 

/(g) = — — / Tr (xu +tuv)dx 
V/(s) ~ —tv — xu 

{It,I(s)} = {-V,tUt +XUx +U) = 

= — / Tr {tvvx + tv[u, v] + xvux + uv)dx 
dl. > 



dt 



2 / ^^ '■""''■ 



In (a.l) the relations 



dl f 1 

{It, I(s)} -^ = - / Tr {tvvx + tv[u, v] + xvUx + -uv)dx (a.l) 



Tr V [u, u] = 

Tr (vVrr) — Tr (v ) 

^ ''' 2dx ^ ' 

Tr (xvUt H — uv) — Tr (xuv) 

^ "" 2 ' 2dx ^ ' 
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hold. In fact, in terms of 7^. one can write 

i k 

= XIX! ^(^* ■-bk){ak - ai){-ix)kil,k 



which imphes 



i k 

Tr [xVxu) 



— — Tr [xuv) = 2Tr {xvu^) + Tr {uv). 



Then: 



dt 2 dx 



{/*,/(,)} ^ = -77 / — Tr {xuv + tv^) dx. 
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